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We use a chiral random matrix model to investigate the 
effects of massive quarks on the distribution of eigenvalues of 
QCD inspired Dirac operators. Kalkreuter's lattice analysis 
of the spectrum of the massive (hermitean) Dirac operator for 
two colors and Wilson fermions is shown to follow from a cubic 
equation in the quenched approximation. The quenched spec- 
trum shows a Mott-transition from a (delocalized) Goldstone 
phase softly broken by the current mass, to a (localized) heavy 
quark phase, with quarks localized over their Compton wave- 
length. Both phases are distinguishable by the quark density 
of states at zero virtuality, with a critical quark mass of the 
order of 100-200 MeV. At the critical point, the quark den- 
sity of states is given by vq(X) ~ l-M 1 / 3 . Using Grassmannian 
techniques, we derive an integral representation for the resol- 
vent of the massive Dirac operator with one-flavor in the un- 
quenched approximation, and show that near zero virtuality, 
the distribution of eigenvalues is quantitatively changed by a 
non-zero quark mass. The generalization of our construction 
to arbitrary flavors is also discussed. Some recommendations 
for lattice simulations are suggested. 
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1. Introduction 



In QCD, the character of the Dirac spectrum near 
zero virtuality follows from symmetries alone, in the limit 
where the quark mass is taken to zero. This is best ex- 
emplified by the Banks-Casher relation |Q for the quark 
density of states, and its moments 0. A number of re- 
cent studies have confirmed these relations using a chiral 
random matrix formulation, and unraveled their generic 
structure around zero virtuality (microscopic or meso- 
scopic limit) in terms of a universal spectral density [§,1). 

Given the importance of the concept of current quark 
masses for spontaneous symmetry breaking and restora- 
tion in QCD, an important yet unanswered question in 
the context of chiral random matrix models, has to do 
with the role of a finite current quark mass. The 1/N 
expansion used in the Coulomb gas approach calls for 
non-trivial subleading effects H[|. The method of or- 
thogonal polynomials Q fails in the massive case || . Re- 
cently, Kazakov 0, Brezin, Hikami and Zee (BHZ) ||, 
and Zee Q have discussed a number of alternative and 
powerful methods to analyze the spectral densities and 
level correlations of a large class of random matrix mod- 



els. Their methods borrow from exact integral represen- 
tations, Grassmannian formulations and diagrammatic 
techniques. Some of these methods will be taken up in 
this paper and applied to various chiral random matrix 
models as inspired by QCD spin and flavor symmetries. 
In section 2, we go over the general aspects of the spec- 
tral density, its relation to the quark condensate and its 
analogy with the Kubo formula. The striking violation of 
Lifshitz's bound JlC) ] near zero virtuality is also discussed. 
In section 3, we recall results from chiral random matrix 
models both in the microscopic and macroscopic limit. 
In section 4, we analyze Kalkreuter's quenched lattice 
SU(2) simulations of the spectral density using recently 
developed arguments by Zee ^ and Kazakov 0. We 
show that for a current mass m = —(N/V^^qq)^ 1 (where 
N/V4 is the number of quark states in the four Euclidean 
volume, and (qq) the quark condensate), the spectral dis- 
tribution for the unsquared and hermitean Dirac opera- 
tor for Wilson fermions, shows a phase transition from a 
Goldstone phase to a phase where the quarks are local- 
ized over their Compton wavelength, with strong chiral 
symmetry breaking. The order parameter in this case 
is the quark density of states at zero virtuality, and the 
transition is reminiscent of a conductor-insulator (Mott) 
transition. In section 5, the unquenched problem is an- 
alyzed using the Grassmaniann method introduced by 
BHZ. A general integral representation for the spectral 
density is derived. In section 6, we show that the integral 
representation leads to the known spectral density in the 
macroscopic as well as microscopic limit for Np = 1 and 
zero quark mass. In section 7, the case of non-zero quark 
mass is investigated, leading to a new spectral density in 
the microscopic limit and for quark masses m N ~ 1 . The 
generalization to an arbitrary number of flavors is out- 
lined. Our conclusions and recommendations are sum- 
marized in section 8. 



2. QCD Spectral Distribution 

The spectral representation for the quark propagator 
of flavor F in a fixed gluon background is given by 



Sf(x, y, A) = ^2 



(1) 



-A„ — imp 

where </>„, A n are eigenvectors and eigenvalues of the 
massless Dirac equation 

if{A)4> n {x) = \ n (t> n (x). (2) 



1 



The fermion condensate in Euclidean space is 

N F 

(Jq)=-J2((TrS F (x,x,A))) 

F=l 



(3) 



where ((. . .)) denotes the averaging over the gluonic con- 
figurations A using the QCD action with massive quarks. 
In the limit where the four- volume V4 goes to infinity, the 
spectrum becomes dense and we may use the eigenvalue 
density §) 



i\,m F ) = {(£5(\-\ n ))). 



(4) 



to characterize the Dirac spectrum. Through the fermion 
determinant in the averaging measure, (Q) carries a non- 
trivial dependence on the current masses m F . We will 
refer to the mass dependence in (|4[) as the sea mass de- 
pendence. In terms of (||) the fermion condensate (||) 
becomes 



\n\ - 



N F 



F=l 



v(\m F ) 
A + im F 



(5) 



The explicit mass dependence in the denominator of (j^) 
will be referred to as the valence mass dependence, for ob- 
vious reasons. In QCD, both the sea and valence masses 
are the same. Here, we may choose to disentangle them 
(for theoretical clarity) whenever indicated. Throughout, 
we will think of the masses as fixed external parameters, 
although in QCD the quantum averaging forces them to 
run. This brings about the nasty issue of the ultra-violet 
sensitivity of (||) for finite current quark masses. Al- 
though perturbative renormalization of is possible we 
will not discuss it here. Most of the discussions to fol- 
low, focuses on the infrared part of the spectrum (around 
zero virtuality). The random matrix models to be dis- 
cussed below are inspired models for the constant quark 
modes only, with some astonishing resemblance to lattice 
regulated simulations. 

In the chiral limit, mp — > 0, the Dirac operator iJJ 
anticommutes with 75, so the non-zero eigenvalues come 
in pairs (A, —A) and the spectral function is symmetric. 
Thus 



irNp 



1/(0). 



(6) 



following a Wick rotation to Minkowski space (q* , q) — > 
(iq, q) . This relation was first derived by Banks and 
Casher It is important that the chiral limit is sam- 
pled with a valence quark mass mp that is taken to zero 
after the thermodynamical limit V4 — > 00, for otherwise 
the result would be zero. The spontaneous breakdown of 
a continuous symmetry cannot take place in finite vol- 
umes, unless the condition m F (qq)V4 » 1 is fulfilled. 
The result (^) states that in vector-like theories with chi- 
ral symmetric (even) spectra, the quark condensate is 



related to the mean spectral density at zero virtuality 
(A = 0) . The delocalization of the quark modes is caused 
by strong correlations that randomize the Dirac spec- 
trum near zero, triggering a huge accumulation at zero. 
As V4 — > 00 the number of eigenvalues grows with the 
four volume V4 in contrast to the length A \/Vl in the free 
case@,0. 

The change in the number of quark states near zero 
virtuality is drastically different from what is expected 
from Lifshitz's condition Jl(],[l2j for the case of scattering 
off random repulsive centers, where the density of states 
is found to vanish exponentially. The reason may be 
traced back to the chirality structure of the random en- 
semble discussed here, hence to the spin of the quarks in 
QCD. Lifshitz's condition can be evaded by noting that 
in a magnetic field further delocalization can be gener- 
ated without cost of energy O]. For spinless particles 
the density of states is bounded from above by the free 
quark density of states Jl| , hence no condensate is allowed 
to form in scalar QCD. 

The Banks-Casher relation (^|) is reminiscent of the 
conductivity in metals, where the latter is proportional 
to the density of states at the Fermi surface. Indeed, it 
follows from the Kubo formula that the d.c. conductiv- 
ity a relates to the density of states at the Fermi level, 
p(Ep), through 



D p(E F 



(7) 



where e is the electron charge and D the diffusion con- 
stant. This result (Q) is reminiscent of (0) with the iden- 
tification of 7rcr/e 2 D with —(qq), and Ep with A ~ 0. 

3. Random Matrix Model 

The pertinent random matrix model for the QCD spec- 
trum near zero virtuality follows from the color represen- 
tation of the quark fields (here fundamental), the chi- 
rality odd character of the massless Dirac operator in a 
fixed gluon background and the number of flavors. In the 
continuum, the QCD Dirac operator for three colors may 
be mapped onto the Gaussian Unitary Ensemble (CUE), 
while for two colors it may be mapped onto the Gaussian 
Othogonal Ensemble (GOE) 0. In the zero topological 
charge sector, and for three colors the generating func- 
tional is §,||J||]i|, 



N F 



Z[m] = dT Yl det 



T 

imp 



-ATE 2 Tr(T t T) 



(8) 



where T is a random N x N complex matrix, with N iden- 
tified with the four volume V4 and £ the chiral condensate 
appearing in the Banks-Casher relation (j^). Equation (j^) 
is the generating function for the GUE. The joint eigen- 
value density following from (||) reads j| 



v(X 1 ,...,X N ;m F ) 
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i<j 
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(g) spectral density. Variation of ( |1 1| ) with respect to v and 
differentiation with respect to A, yield 



with an overall normalization (see below). Inte- 
grating (||) over (N — l)-eigenvalues and taking the mi- 
croscopic limit N — > co with x = NX hxed||] yields 
the following form for the microscopic spectral density 

K = o) [§ 

v s {x) = 2Y?x( J 2 (2Ex) - J A r F _ 1 (2Sa ; )JAr F+1 (2Ex) ) (10) 



which is the master formula for all the sum rules dis- 
cussed by Leutwyler and Smilga [0 . Equation (|H]) shows 
that around zero virtuality (A ~ 0), the distribution of 
eigenvalues oscillates to zero. For Np = these oscilla- 
tions are caused by the level repulsion around zero due to 
the symmetric character of the spectrum under chirality 
(Airy phenomenon ||). These oscillations are affected 
by the fermionic determinant in the massless case, as is 
clear from the Np dependence (zero mode suppression). 
In the quenched approximation, the character of these 
oscillations appear to be unaffected by the choice of the 
measure, provided that it is local. Non-local changes to 
the measure, e.g. through a fermion determinant, do af- 
fect the structure of these oscillations. Other non-local 
changes are also possible, but will not be discussed in 
this work. 

The joint eigenvalue density associated to (0) can be 
obtained through a direct integration, or by using a 
Coulomb gas description of (||). Indeed, using the def- 
inition (Q) for v(X), where the average is over the joint 
eigenvalue density (^|), allows a rewriting of (|J) in terms 
of an effective action 

S[u] = - J dXdX 1 u(X)v(X') ln|A 2 - A' 2 | 

- J dv v(X) ( ^ ln(A 2 + m%) + ln|AM 
+NY 2 [ v(X)X 2 



(11) 



+a / dXv(X)-N 



where £ is a Lagrange multiplier. For a dense spectrum, 
the integration over the eigenvalues A; may be traded 
by a functional integration over the eigenvalue density 
v{X), modulo a Jacobian ||. In the large N limit, the 
extremum of S[^] including the contribution from the Ja- 
cobian to the effective action, determines the macroscopic 



2P 



u{X') 



TVS 



y L_ 

z -' A 2 + m 

F=l 



A 2 - A' 2 

with the normalization condition 
dX v{X) = N 



1 

2A 2 



-J (12) 



(13) 



The contribution J is due to the Jacobian and is of order 
IN 2 ||. Its explicit form will not be needed here. In 
P stands for the principal value of the integral. In 
the thermodynamical limit [N — > oo), the fermions and 
the Jacobian drop from the macroscopic spectral density 
in the chiral limit. Thus 



S 2 

= N — 
A 2 - A' 2 2 



iz(A') 



(14) 



This integral equation yields a semi-circular distribution 
for the macroscopic spectral density 



NY, A 2 £ 2 

v{x) = — v 1 - — 



(15) 



The level repulsion revealed in the microscopic limit 
takes place in a window of size 1/N around the ori- 
gin, and shrinks to zero size in the thermodynamical 
limit. In terms of (|l5|), the quark condensate is (qq) = 
—Np(N/V4)T,. Here N/Va is just the eigenvalue density 
of the Dirac operator in Euclidean space. Throughout £ 
and V4 are set to one, and the thermodynamical limit is 
understood for N — * 00. The scale E can be reinstated at 
the end by inspection. Since the thermodynamical limits 
are now understood, we prefer to work from this moment 
with the macroscopic spectral density ( |l5| ) normalized to 
one instead of N, i.e., 



u(X) 

unless specified otherwise 



(16) 



4. Quenched Spectral Distributions 

Recent detailed numerical analysis by Kalkreuter fltt l 
using Wilson as well as staggered fermions, provides 
some useful insights to the macroscopic character of the 
fermionic spectrum of four-dimensional gauge theories. 
On the lattice, the spin-Lorentz structure of the QCD 
Dirac operator is affected. Indeed, for any number of 
colors greater than twoQ and Wilson fermions the per- 
tinent random matrix ensemble is the Gaussian Unitary 



1 Note that our definition differs slightly from the one used in 
reference Q due to replacement 2N <-> TV, therefore v s (x) 
2u s (2x). 



2 For two colors the ensemble is the Gaussian Orthogonal 
Ensemble (GOE). 
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Ensemble (GUE). The chiral structure is upset by Wil- 
son's r-terms, needed to remove the lattice doublers. For 
two colors and staggered fermions the chirality structure 
is preserved, but the Lorentz structure is upset. The 
pertinent random matrix ensemble is the Gaussian Sym- 
plectic Ensemble (GSE) ||. For three colors it is back to 
the GUE. 

• Wilson Fermions 

In the quenched approximation, Kalkreuter's results 
for SU (2) c with Wilson fermions for the unsquared and 
unnormalized operator Q = 75 (JJ + m) (continuum), 
may be mocked up by the Gaussian Orthogonal Ensem- 
ble (GOE), provided that the gauge configurations are 
sufficiently random. Specifically PL 



Qw 



m 

— m 



R 



(17) 



which is the sum of a deterministic 2N x 2N matrix (first 
contribution, with m - diagonal block N x N), and a 
random hermitean (symmetric) 2N x 2N matrix. We 
note that QL, = Qw- The measure is 



P(R) = — e -2ATTr(RRt) 

z 



(18) 



This problem is reminiscent of an electron scattering on 
impurities in a spin dependent quantum Hall fluid, as 
recently suggested by Zee (9). 

The distribution of eigenvalues of Qvy 

v Q (X,m) = ^ <Tr 2N 8(X-Q w )> (19) 

follows from the resolvent G(z,m) of Qw through 



vn(X 1 m) = Im G(A + iO, m) 

7T 



(20) 



Since R is only hermitean (symmetric), the distribution 
of eigenvalues is only symmetric about zero virtuality 
on the average. In the chiral limit, the states are not 
necessarily paired about zero, because of the r-terms in 
the Wilson action. 

The problem of determining the one-point Green func- 
tion G(z) for the sum of a deterministic Hamiltonian Hq 
with eigenvalues £i, (i = 1, . . . , N) and a Gaussian ran- 
dom matrix was first solved by Pastur [jl7|, and recently 
rederived and generalized using much simpler arguments 
by Zee Generally, the Green function is determined 
from 



G(z) 



1 N 



1 



G(zY 



(21) 



For the case (|17|), the deterministic hamiltonian is com- 
posed of two diagonal blocks with eigenvalues ±m, so 
G(z,m) satisfies the cubic equation 



G 3 - 2zG 2 + (z 2 - to 2 + 1)G - z = 
Using Cardano's complex solution to ( p2| ) we get 

v Q (\m) = ^ [(r + Vd) 1/3 - (r - 
with 



r = -X ( 1 + 2to 2 - -A 2 
6 I 9 



(22) 



(23) 



d = 



1 

27 



(1 - m 2 ) 3 - A 2 (- - 5m 2 - 2m 4 + to 2 A 2 ) 



(24) 



For one flavor we set tuf = m. 



At the critical point ~ 1 and near zero virtuality A ~ 
0, the order parameter behaves as vq(X, 1) ~ IAI 1 / 3 , with 
the critical exponent (3 = 1/3. This is to be contrasted 
with the density of states of electrons near the mobility 
edge, rescattering off random impurities and undergoing 
Anderson localization with j3 — Jl8|-p0[. 

It is worth mentioning, that although the random ma- 
trix model under consideration is only valid for the con- 
stant quark modes, it does reproduce the bulk character- 
istic of the spectral function from lattice regulated cal- 
culations. This justifies a posteriori our assumption in 
ignoring the ultraviolet aspects of the quark spectrum, 
with their inherent diverging contribution to the quark 
condensate. We suspect that a cooling of Kalkreuter's 
gauge configurations will not affect considerably the char- 
acter of the spectral distribution. Such procedure can be 
used to define unambiguously the quark condensate for 
finite current masses. 

The distribution of eigenvalues of Qiy is strikingly sim- 
ilar to the one discussed recently in the finite tempera- 
ture problem for the hermitean and massless Dirac op- 
erator ip on the torus in the quenched approximation 
[f2~l| [23I . This may be understood if we note that in the 
high temperature phase, dimensional reduction implies 
that in Euclidean space and in one-dimension lower, the 
effective quark mass is not m but V m 2 + ir 2 T 2 after suit- 
able chiral rotations [p4| . In three dimensions the ran- 
dom ensemble is indeed the unitary ensemble p5] . This 
observation can be used to map out the chiral sensitive 
part of the QCD phase diagram with Wilson fermions 
in the plane spanned by the temperature and current 
quark mass. We note that in the context of phase tran- 
sitions, and if G (or at least its discontinuity along the 
real axis) is understood as an order parameter, then the 
cubic equation (|2^) is generic of second order phase tran- 
sitions, as expected from universality arguments. Some 
of these issues will be discussed elsewhere. 

Kalkreuter's results are displayed in Fig. [I] for (3 — 
A/g 2 = (strong coupling) and different values of n = 
(2m + 8) _1 (Wilson fermions). A similar, although 



4 



weaker, behaviour is also seen for j3 = 1.8 (weak cou- 
pling) and for a range of k's close to the critical value of 
0.125 M. 





FIG. 2. Spectral function ttvq(\,iti) (|23|) as a function of 
the mass m and eigenvalue A. Point A fixes the normalization 
for the massless case, point B corresponds to the critical mass 
m, — 1. All axes are scaled in units of E = 1. 





FIG. 1. Spectral densities for the operator Qw for Wilson 
fermions at /3 = 0, obtained by Kalkreuter [16]. The stars 
(upper), crosses (middle), squares (lower) correspond to the 
values of k equal to 0.25 ,0.20, 0.125, respectively. 



The behavior of ( |23| ) versus A for different values to is 
shown in Fig. @. 



In the massless case, ( j23| ) reduces to ( |l5| ) in units of 
1/S, with a half-width equal to 2. In the limit m >> 2, 
the density of states decompose into two pieces, with 
i/q(0, to) = 0, and a width of the order of 3/^/2 which 
is to be contrasted with 4, the size of the semi-circle. 
This is a regime, where chiral symmetry is strongly bro- 
ken by a massive quark. For to* = 1, the two regions 
merge into each other. In physical units, the spectral 
transition occurs for quark masses of the order of to* ~ 
-(N/Vi) {uu}' 1 . Typically N/Vi ~ 1 fm" 4 , and for a 
quark condensate in the range (200MeV r ) 3 - (250MeV) 3 
we get a critical mass in the range to* <~ 100 — 200 
MeV, which is rather close to the strange quark mass. 
The spectral transition characterizes the transition from 
a delocalized phase with softly broken chiral symmetry 
(Goldstone phase) to a localized phase with strongly bro- 
ken chiral symmetry (heavy quark phase). The transition 
is characterized by the quark density of states at zero 
virtuality, ~ 0,to). It is non-zero in the Goldstone 
phase, and zero in the heavy-quark phase. This transi- 
tion is reminiscent of a Mott-transition from a conductor 
(Goldstone phase) to an insulator (heavy quark phase). 

• Staggered Fermions 

For staggered fermions and two colors, it is more 
appropriate to use the Gaussian Symplectic Ensemble 
(GSE) H|, while for three colors the GUE. For illustra- 
tion, consider the GUE for three colors, that is 



TO 

— m 



-iT 
iTt 



(25) 



which is again the sum of a deterministic 2N x 2N ma- 
trix (first contribution, with to - diagonal block N x N) , 
and a random hermitean complex 2N x 2A^ matrix (sec- 
ond contribution, with T - block N x N). The resolvent 
associated to (E5J) is 
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. z I V4 + m 2 -z 2 

G-(z) = -11 — % — ===== — 

2 \ Vz 2 — to 2 



(26) 



and the corresponding macroscopic spectral distribution 
is 



vq (A, to) = 



|Aj V4 + m 2 -A 2 

2-7T VA 2 — TO 2 

x9(A 2 -to 2 ) 9(4 + to 2 -A 2 ) (27) 



where 9 are step functions. The macroscopic spectral 
distribution is symmetric about the origin, with a sup- 
port from m to y4 — m 2 to the right, and — v 4 — m 2 and 
— m to the left. The appearance of the gap at zero vir- 
tually is remiscent of the familiar mass gap in the Dirac 
equation. Indeed, Qs is nothingh but the Hamiltonian 
of a quark in a five-dimensional space where 75 plays the 
role of the fifth gamma matrix along the extra (temporal) 
direction. 

The spectral distribution for the squared staggered op- 
erator Q|, 



Q 



2 / m 2 + TTt 



m 2 + T+T 



that is 



v Q2 (\ 2 ,m) = ^< Tr 2W (A 2 - Q|) > 



(28) 



(29) 



with the averaging carried out using ( |18| ) , can be readily 
tied to (E7D through 



1 



i/ Q3 (A ,m) = — v Q (X,m) 



(30) 



This result can also be checked to follow from the large 
N analysis using the method suggested by Kazakov or 
BHZ. 

To show this, consider the modified density of states 



z>(A, m) = VA 2 - to 2 « Tr N 5{\ 2 - m 2 - T f T) >> (31) 

which is seen to reduce to (|l^) in the massless case, 
£(A, 0) = v{X). Note that the normalization is now to 
N. In terms of (|3l|), the quark condensate reads 



m 1 



IaIVa 2 " 



9(X,m)dX (32) 



Following Kazakov [|7| , we modify the Gaussian probabil- 
ity distribution (|l^) by adding an auxiliary matrix source 
A, 



Pjt(T) = 1 c-NTrjT^T-AT^T) 

Za 



(33) 



As a result, the Fourier transform of the resolvent reads 
U a {t) = VA 2 - to 2 ^l T re lrTtT+trm2 ^ (34) 



where the subscript A denotes averaging with the mod- 
ified probability distribution (|33|). In terms of ([m]), the 
spectral density is 



to) = y/\ 2 - m 2 I ^e trX ' 

, 27T 



£/a=o(t) (35) 



where A is set to zero, only after the averaging in ( |34| ) 
is carried out. The result for the macroscopic spectral 
density is 



9(X,m) = — \/ 4 + m 2 — X 2 

Z7T 

x9(A 2 -m 2 ) 9(4 + m 2 - A 2 ) (36) 

in agreement with (|30| - ^ll) . 

The above arguments are subtle at the edge of the spec- 
tral distributions where Airy oscillations are expected [0 . 
These oscillations are at the origin of (|l(J) in the quenched 
approximation. To analyze them in the staggered case, 
we need to magnify the edge points at the level of one 
level spacing. This can be achieved for instance, by tak- 
ing the microscopic limit N — > 00 with Ny/X 2 — m 2 fixed. 
For the chiral unitary ensemble, the corresponding micro- 
scopic spectral density reads (quenched approximation) 



N 



P s (X,m) = yvA 2 



x (j 2 (2iV V / A 2 - m 2 ) + J?(2Ny/>? 



(37) 



In the region ±m, the spectral density dives to zero again, 
at a rate that is comparable to the one discussed in the 
massless case (|l(i|). The effects of the mass on the un- 
quenched spectral density near zero virtuality are more 
difficult to track down with the methods described above. 
They may be analyzed by Grassmannian techniques as 
we now discuss. 

5. Integral Representation : Np = 1 

In a recent analysis, BHZ have put forward an alter- 
native method to the Coulomb gas approach and the 
orthogonal polynomial construction to study the non- 
Gaussian character of the density of eigenvalues. Their 
method uses Grassmannian techniques, and proved to 
be very elegant to get at the oscillating character of the 
density of states near zero virtuality. Here, we will use 
their method to analyze the elusive case of massive, un- 
quenched QCD with three colors. Using the GUE, we 
will discuss the one flavor case in details, and outline the 
generalization to the case of two and more flavors. 

The density of eigenvalues i/(A, to) will be again sought 
in terms of the discontinuity of the resolvent G(z, to) 
along the real axis. The latter is now given by 



Tr 



2N 



1 



M 



(38) 
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with 



M = 



T 
Tt 



(39) 



The averaging in ( p8[ ) is over the distribution P(T) for 
complex N x N matrices T 



P(T) 



Z 



N F 

] [ det 2N M F 



(40) 



F=l 



where Z is an overall normalization, including the 
fcrmion determinant. For each flavor, 



M, 



im T 
im 



(41) 



First, let us consider the case of only one flavor Np = 1. 
Following the method discussed by BHZ, we can write the 
unaveraged contribution to the resolvent as follows 



det 2 7v Mi 
= J dn({a\a) + {b\b)) 

a iNz((a\a) + (b\b)+( a \a) + (f3\P)) 



xe 



xe 



-iAf«a|T t |()) + <Q|T t |/3)-< At |T t |i/>+h.c. 
-Nm{{^) + {v\v)) 



(42) 



where di,b{ are N-dimensional complex vectors, and 
Q!i, are N-dimensional Grassmannian vectors. 

The measure dfi is graded, 



N 



dn = Y]_ [da-i] [dbi] [den] [d/3i] [dfj,i] [dui 



(43) 



Averaging (42) over T using the Gaussian measure ( ]40| ) 
gives, 



1 r N 

G(z,m) = - J Y[dfi((a\a) + (b\b)) 



JNz ({a\a) + {b\b)+{ a \a) + {f3\P)) 
v<) -JVto((a 1 |m>+Wi'»--W'&n 



with 



Tr JV = Tr JV ((-|6){a| + | ( 9)(a|+H(//|) 
®(-|a)<&| + |a)</3| + |A*>M)) 



(44) 



(45) 



The tensor product in ([44] ) generates four Fermi fields 
bilinears. They can be linearized by introducing four 
auxiliary complex fields <Jij with i,j — 1,2, e.g. 



e N(<x\a)W) = 

N f J2 

H / d erne 

7T ./ 



A r (<T 1 * 1 <Tii+<T 1 * 1 <a|a)+o- 11 </3|/3)) 



(46) 



and so on. In terms of the Grassmannian integra- 
tions can be undone, and the result is 



N 



G(*) oc / JJ[dai][d6i][do-]((o|o) + (b\b)) 



1=1 



iJVz((a|a>+(6|6»-iV<a|a>(6|6> 



xe 
xdctA 

where AN x 4A^ block matrix reads 

/ *i 

A= \b)(a\ 
\ 



a)(b\ 


• s 3 





Sl 





s 4 





So 

\b)(a\ 


|a)(6| 


«3 


S2 



(47) 



(48) 



with s\ — iz — an, s 2 — ra — 022, S3 = — <j\i and S4 = 
— (72i. The matrix ( |4S| ) is sparse. It can be rotated to an 
block-diagonal form of 4 x 4 matrices. Hence 



detA = {{a 2 b 2 f - a 2 b 2 ^ Si s* + \ Sl s 2 - s 3 s 4 \ 2 )) 



i=l 



X [\S1_S2 ~ S 3 S 4 \ 



,N-1 



(49) 



with a 2 — (a\a ) and b 2 = (b\b). The integrals in @ may 
be simplified by introducing the new variables 



izSn 
izSn 



m8 2 2 ■ 



(50) 
(51) 



and x = a 2 and y = b 2 . Sn is non-zero for i = j = 1 and 
similarly for S 22 . Thanks to the Gaussian term in the 
integration over a and z = A + ie, we can freely shift the 
contour of integration, in particular we can treat matrices 
s and s* as complex conjugate. Hence 

/>oo 

G(z, m) = Af e - Nm2 / dx dy ds tj ds* tj 
Jo 

X(X + y)( X y) N ~le- NX V- NtlSS ' +Nm{s ^ +S ^ 

x e iC(x+J/-8n-^i) e C 2 /w detA(ss t ,a;y) (52) 
where Af is a normalization constant and 



det A = (xV 



xy trs s T + det ss') (det s s T ) 



(53) 



Notice that the saddle point in Sij has a symmetry with 
respect to the unitary rotations of s for m = 0. Only the 
O(l) term breaks this symmetry. The relations (52 5^) 
constitute the integral representation for the resolvent 
G(z,m) in the complex z-planc, for one massive flavor 
and finite dimension N. Its singularity along the real axis 
z = A + ie corresponds to the distribution of eigenvalues 
of the QCD Dirac operator with one massive flavor as 
modeled by the chiral random matrix ensemble. 

In the thcrmodynamical limit N — > 00, the integral 
in ( p2| ) is dominated by the saddle point configurations. 
The latter are the extrema of 
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5cff = xy - In xy - iz(x + y) 
+Tr(ss t ) - ln(detss t ) 



(54) 



The extremum of ( p4[ ) coincides with the extremum in 
the quenched approximation (BHZ), leading to a semi- 
circular distribution of eigenvalues. The fermion deter- 
minant does not affect the eigenvalue distributions in the 
large N limit. As observed in j^] and discussed by BHZ, 
the saddle point approximation breaks down at A = 
and also the edges of the semi-circular distribution due to 
Airy-type singularities (quenched) and zero modes (un- 
quenched). 

6. Microscopic Limit : Massless Case 

To investigate the behavior of the spectral density near 
the origin for the GUE, we will use the microscopic limit 
discussed in ||. We will simply blow up the spectrum 
near A = over a scale of the order of one-level spacing. 
This is achieved by considering Nz — ( = 0(1) in the 
region z = (zero virtuality), and track down the l/N 
corrections. Since the mass term and the l/N correc- 
tions violate the saddle point invariance under unitary 
rotations, care is to be used to disentangle the soft direc- 
tions from the hard directions. All soft directions have 
to be integrated out exactly, for an accurate estimate of 
the saddle around z = 0. 

First, let us consider the massless case m = 0, where 
only the 0(1) terms break the unitary symmetry of s. To 
make the latter manifest, we use the parameterization as 



S21 



so that 



8\2 
.S'22 



HI 







Jrf> 




i/3 



cost/ 
— sin( 



sin# 
cos 9 



trss' = u + v + w, 
det ss' = uw 
s n + s n = 2v / ucos#cos< 

Following BHZ let us introduce also 

xy = p, 
x + y = 2 v /pcoshw. 



(55) 



(56) 
(57) 
(58) 



(59) 
(60) 



With these notations we have neglecting the 0(l/N) 
term in the exponent 



/*°° /*°° Ar\ 

G(z,0)=N J / coshwdw / 

Jo Jo VP 



/>00 p27T 

x / ududvdw / d(f)di[)dl3d'y (61) 
Jo Jo 

,tt/2 

x / sin 9 cos 9d9p^ {uvS)^ 



x (p 2 — p(u + w + v) + uw) 

x g— N (u+w+v+p) 

X g^^CVP c °sh uj^2iQ^/u cos (j> cos 9 



(62) 



(63) 



The tp,(3 and 7 angular integrals can be trivially per- 
formed. Also the v integral is trivial. The 9 and <fi inte- 
grals can be performed using 



2ir i-it/2 



d9 sm9 cos9e- i<: ^ lcos9cos ' l> 



Jo 
Ji(2Cy^) 

and the u integral by 



/ duo coshwe 2iC ^ cosh " = J ftT 1 (-2<^) 
Jo 



(64) 



(65) 



;(Ji(2C y /p)+iN 1 (2C y /p)). 



For the imaginary part of G(£, 0) we get 



e"" l — du — 

Vp w 

1 



X {P — P( U + W + Jr) + UW ) 

x (puw) N e- N ^ + ^ M2CVP) ■ ( 66 ) 

Where we have explicitly wrote down the normalization 
(we used Stirling formula to rewrite the normalization in 
this form). This integral is dominated by the neighbor- 
hood of the saddle point 



Writing 



p c = U c = W c = 1. 



P= 1 + ^ 

VN 



u=l + 



w = 1 + 



(67) 



(68) 



N 



N 



we have 



p 2 — p(u + w + 1 /N) + uw — 

l((p') 2 - 1 - u'pi + «/(«' - p>) - -Lp>). (69) 

The leading contribution coming from (p') 2 — 1 cancels 
out. We are therefore led to consider the next order 
terms. We observe that due to the symmetry in p' and 
u' the terms containing w' cancel. Using (K39) we have 
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Ji(2Cv^) 



Let us introduce 



u 1 (u')\ 2 



t = cos ( 



^(20-^=^2(20 



Fn 



2N 



u' {u') 2 



(70) 



log u ) 



-N(u— log u) 



Using this expansion we get 

p((/N) = -2CJi(2C)J 3 (2C) + Ji(20 J 2 (20 + CJ 2 2 (20 
= 2C (^(20-^(20^(20)- (72) 

in agreement with (p~0|) for Np — 1. 



7. Microscopic Limit : Massive Case 

Let us consider the general case of one massive flavor. 
The inclusion of the mass corresponds to analyzing (52). 
Three different regimes may be be considered, depend- 
ing on the strength of the sea quark mass. The regime 
Nm » 1, corresponds to the case where the sea mass 
is comparable to the valence mass, in thephysical (chi- 
ral) limit. In this case, the mass term in ( p2[ ) contributes 
to the saddle point equations along s and s'. The sad- 
dle point equations in x and y remain unaffected, and 
the resulting spectral distribution in the large N limit 
is a semi-circle. This result is in agreement with the 
Coulomb gas analysis. We expect, the microscopic spec- 
tral density to be unaffected by the mass effect. In the 
regime Nm << 1, the mass effects are sub- leading com- 
pared to the 0(1) terms. This regime corresponds to 
the massless case discussed above. Both the macroscopic 
and microscopic spectral densities are expected to be un- 
affected by such masses. The most interesting regime, is 
the one for which Nm ~ 1, as we now discuss. 

In the case ji — Nm = 0(1), 



G(z,m)=J\f j dx dy dsij ds*j(x + y)(xy) (73) 



x e 



^ 2 / N det A(ss t ,a; 2 /). 



We can use the same parameterization as before and ob- 
serve that 

S22 + s* 22 = 1\fv sin 9 cos(/3 + ip) + 2\fw cos 8 cos(7 + ip). 

(74) 



We have 

G{z,m)=M > 



r<>s!i„vL: / -1— 

o Vp 



(75) 



(76) 



X / ududvdw / dcftdipdpd-/ 
let Jo 

„JV t„.„..\N-l/Jl 



2tt 



x / t dtp" (uw) N 1 (p z -p(u + w + v) + uw) 
o 

- N (u+w-\-v+p) J2i<^ ^/p cosh uj 2iC,^/ut cos (p 



X e 



— [i^/wt cos 7— fly/v^/l — t 2 COS /3 



We perform the integrals over u> (as before) and over 
angular variables (j), ip, (3 and 7 using 





2tt 



e lxcos4 > = 2vrJ (x), 
d<j) e' xcos ^ = 27r/ (aO- 



(77) 
(78) 



We get 



p((/N,fi/N) =-n 2 V2 



N 



7/2 



^3N 



dp dv du 



dw 



t dt(p 2 — p(u + w + v) + uw) 



(79) 



x (puw) N e~ N ^+ u+w+ ^ Jl(2C » j (2Ctv^) 

Vp 

x 7 (2/iVl - t 2 Vv)h{2^itVw). 

This is an explicit integral representation of the spectral 
distribution in the massive case and for finite N. No ap- 
proximation has been used in getting from (f74|) to (j8fj). 
We observe that the integral over v can be performed 
in the leading order and that to this order we can re- 
place Io(2/i\/l — t 2 y/v) by one. Using this fact and the 
formula: 



1 1 dtJ (xt)I (yt) = vKv)M4 +xM*)Io{v) m 



x 2 - + y l 



we have 



/•oo dw 1 

p(0^)oc / dp du — (p 2 — p(u + w + — ) + uw) (81) 
Jo w N 



(puw) 



N —N(p+u+w) 



Vp 



At v /wJi(2/xy^) J (2Cv^I) + CV^Ji(2Cy/n)I {2^) 
X 2((0 2 m + p?w) 

The integral is dominated again by the neighborhood of 
the saddle point at p c = u c = w c = 1. We need as before 
terms 0(1/ N) to calculate it. 







Expanding the above formula as before, we obtain after 
somewhat lengthy algebra, 



3 + cos 4t 



p(C,n) = 2I Jq A cost - J Ji 



- J 



-8A 2 cos t + cos 3t — cos 5t 



SA 



2I\ sin r J cos t + JqJ\- 



cos2rsin r 



A 



2 J? cost sin 2 t) 



with the parameters 



C = A cos t, 
/.j = Asinr. 



(82) 



(83) 



The Bessel functions in ( p2[ ) are defined as J„ = J„(2£) 
and /„ = 7 n (2/i). For r = 0, we recover ( |ic| ) for JVj? = 1 
and 

Clearly, sea masses of the order Am = 0(1) do affect 
the character of the spectral oscillations near zero virtu- 
ally. This in turn would imply, new sum rules for the 
moment of the eigenvalues of the massive QCD Dirac 
operator. The very non-local character of the fermion 
determinant upsets the universality at zero virtuality. 

The present analysis can be extended to arbitrary fla- 
vors Np = 2,3, .... Indeed, for N F flavors the bosoniza- 
tion of the fermion bilinears requires the introduction of 
2(Np + l) 2 auxiliary fields, which we label generically by 
the matrices s and . After integration, the analog of the 
determinant A in (§§) splits into (2N F + 2) x (2N F + 2) 
blocks each of size N x N. By analogy with the one fla- 
vor case, the determinant can be reduced to the product 
of N determinants of size (2N F + 2) x (2N F + 2). The 
result, can be generically written in the form 



'TV* 



P A r F+1 ((a 2 6 2 ),sst)det Ar - 1 (sst) 



(84) 



where P is a polynomial of degree N F + 1 in a 2 b 2 with 
coefficients depending on the invariants of the combina- 
tion ss^ . So the number of flavors enters explicitly as the 
degree of the polynomial and implicitly as the dimension 
of the matrix s. In the large N limit the integral is dom- 
inated by the saddle point, in which the polynomial does 
not contribute. Hence, the saddle point configurations 
are the minima of the action (54), resulting into a semi- 



circular distribution of eigenvalues. This result is totally 
in agreement with the Coulomb gas argument. We ex- 
pect, the above arguments near zero virtuality to carry 
(tediously) to the massive case. 

8. Conclusions 

Using chiral random matrix models as inspired by 
QCD spin and flavor symmetries, we have discussed the 
effects of a light quark mass on the QCD Dirac spectrum. 
In the thermodynamical limit (large N), the quarks play 



a subleading role on the macroscopic spectral density, ex- 
cept near zero virtuality. For one flavor and three colors 
(GUE), we have shown that the microscopic spectral dis- 
tribution is affected by quark masses Am = 0(1). The 
non-local character of the fermion determinant causes the 
spectral oscillations to be flavor and mass dependent near 
zero virtuality. New sum rules of the type discussed by 
Leutwyler and Smilga are therefore expected. Our argu- 
ments extend to more than one flavor. 

Kalkreuter's spectral distributions from lattice SU(2) 
gauge theories make use of a hermitean variant of the 
Dirac operator. In the quenched approximation and for 
Wilson fcrmions, they can be understood using a Gaus- 
sian Orthogonal matrix model, reminiscent of the one 
used in quantum Hall fluids. In contrast to the preced- 
ing discussion, here the quark mass survives the ther- 
modynamical limit. For one flavor, Kalkreuter's spectral 
distributions show that a spectral transition occurs for 
quarks masses of the order of the strange quark mass. 
The transition is reminiscent of the Mott-transition from 
metals to insulators, with the density of states at zero 
virtuality playing the role of an order parameter. 

The similarity between the spectral distribution gener- 
ated from lattice simulations and the one obtained from 
the chiral random matrix model with massive quarks sug- 
gests that the lattice results may be robust to cooling. If 
that is the case, then we suggest that a properly regu- 
lated quark condensate in the presence of finite current 
quark masses can be obtained from cooled lattice config- 
urations, in agreement with the general lore of random 
matrix theory. This point is worth checking on the lat- 
tice. 

The behavior of the spectral density near zero virtual- 
ity and for massive quarks, is of relevance to finite tem- 
perature lattice simulations of QCD p^ . The sensitiv- 
ity of the lattice results both to the masses and finite 
volumes (whether temporal or spatial) in the infrared 
regime, are important for a reliable assessment of the chi- 
ral condensate, and hence the study of the spontaneous 
breaking of chiral symmetry in finite systems both in vac- 
uum and matter. The role of the current quark masses 
and the number of flavors in the QCD phase transition 
is of paramount importance for issues such as the order 
and character of the transition. Some of these issues will 
be discussed elsewhere. 



Acknowledgements 

We are grateful to Dr Thomas Kalkreuter for send- 
ing us his data. M.A.N, wishes to thank Gabor Papp, 
for discussions and help with computer graphics. This 
work was supported in part by the US DOE grant DE- 
FG-88ER40388 and by the Polish Government Project 
(KBN) grant 2P03B19609. 



10 



[1] T. Banks and A. Casher, Nucl. Phys. B169 (1986) 103. 
[2] H. Leutwyler and A. Smilga, Phys. Rev. D46 (1992) 
5607. 

[3] J.J.M. Verbaarschot and I. Zahed, Phys. Rev. Lett. 70 
(1993) 3852. 

[4] J.J.M. Verbaarschot, in Continuous Advances in QCD, 
Ed. A.V. Smilga, World Scientific (1994); and references 
therein. 

[5] M. L. Mehta, Random Matrices (Academic Press, N. Y. 
1991). 

[6] A. A. Migdal, Mod. Phys. Lett. A8 (1993) 139. 

[7] V. A. Kazakov, Nucl. Phys. 354 (1991) 614. 

[8] E. Brezin, S. Hikami and A. Zee, "Oscillating density of 
states near zero energy for matrices made of blocks with 
possible application to the random flux problem" , eprint 



cond-mat/9511104. 



[9] A. Ze e, "Law of additio n in random matrix theory" 
eprint |cond-mat/9602146 . 



[10] M. Lifshitz, Sov. Phys. Usp. 7 (1965) 549. 
[11] A. Casher and H. Neuberger, Phys. Lett. B139 (1984) 
67. 

[12] H. Neuberger, Phys. Lett. B112 (1982) 341. 

[13] M.A. Nowak, J.J.M. Verbaarschot and I. Zahed, Nucl. 

Phys. B324, 1 (1989). 
[14] Y.A. Simonov, Phys. Rev. D43 (1991) 3534. 
[15] E.V. Shuryak and J.J.M. Verbaarschot, Nucl. Phys. 

A560 (1993) 306. 
[16] T. Kalkreuter, "Numerical analysis of the spectrum of the 

Dirac operator in four- dimensional 5*17(2) gauge fields", 

eprint |hep-lat/9511009 . 
[17] L. A. Pastur, Theor. Math. Phys. (USSR) 10 (1972) 67. 
[18] F.J. Wegner, Z. Phys. B35 (1979) 208. 
[19] E. Abrahams, P.W. Anderson, D.C. Licciardello, and 

T.V. Ramakrishnan, Phys. Rev. Lett. 42 (1979) 673. 
[20] A.J. McKane and M. Stone, Ann. Phys. 131 (1981) 36. 
[21] M. A. Stephanov, "Chiral symmetry at finite T, the 

phase of the Poly akov loop and the spectrum of the Dirac 

operator", eprint trep-lat/960100l| . 
[22] A. D. Jackson and J.J.M. Verbaarschot, "A random ma- 
trix model fo r chiral symmetr y breaking", SUNY-NTG- 

95/26, eprint |hep-ph/9509324j . 
[23] A. D. Jackson, M. V. Sener and J.J.M. Verbaarschot, 

"Univ ersality near zer o virtuality", SUNY-NTG-96/1, 



eprint hep-ph/9602225 



[24] T.H. Hansson and I. Zahed, Nucl. Phys. B374 (1992) 
117; G.E. Brown et al., Phys. Rev. D45 (1992) 3169; 
T.H. Hansson, M. Sporre and I. Zahed, Nucl. Phys. B427 
(1994) 545. 

[25] J.J.M. Verbaarschot and I. Zahed, Phys. Rev. Lett. 73 
(1994) 2288. 

[26] M. Fukugita, T. Kaneko and A. Ukawa, Phys. Lett. B130 
(1983) 199. 

[27] see e.g. K. Kan aya, "Finite Tem perature QCD on the 
Lattice" , eprint |hep-lat/951004C . 



11 



